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Abstract 

The search for particle electric dipole moments represents a most 
promising way to search for physics beyond the standard model. A num- 
ber of groups are planning a new generation of experiments using stored 
gases of various kinds. In order to achieve the target sensitivities it will be 
necessary to deal with the systematic error resulting from the interaction 
of the well-known v x E field with magnetic field gradients (often re- 
ferred to as the geometric phase effect |1U|. This interaction produces 
a frequency shift linear in the electric field, mimicking an edm. In this 
work we introduce an analytic model for the correlation function which 
determines the behavior of the frequency shift and show in detail 
how it depends on the operating conditions of the experiment. We also 
propose a method to directly measure ths correlation function under the 
exact conditions of a given experiment. 
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1 Introduction 

The proposition that the search for particle electric dipole moments (edm) rep- 
resents a reasonable method to look for physics beyond the standard model 
is inspiring many groups to search for edm's in a variety of systems. (See j5] for 
a recent review). Experiments on several systems including the neutron [3], and 
several species of confined gases ^ including Radium [S] , Radon |B] and Xenon 
\7\ are in various stages of preparation. These experiments are all hoping to 
reach sensitivities in the range of 1CP 27 — 10 _28 e — cm. Sensitivity in this range 
has already been achieved in the case of Hg The experiments proposed 
represent a broad range of operating conditions, from room temperature gases 
with buffer gas to laser cooled atoms in a MOT. 

In order to achieve the target sensitivities it will be necessary to deal with 
the systematic error resulting from the interaction of the well-known 1? x E 
field with magnetic field gradients. Often referred to as the geometric phase 
effect [5], ^U] this interaction produces a frequency shift linear in the electric 
field, mimicking an edm. This systematic effect is highly dependent on the 
operating conditions of the experiment. While experiments in small vessels and 
with high pressure buffer gas are expected to be relatively insensitive to the 
systematic effect, each of the proposed experiments will have to be analyzed 
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in detail to judge its sensitivity to the effect and to find methods of dealing 
with it. In this work we introduce an analytic form of the correlation function 
which determines the behavior of the frequency shift [IT] and show in detail 
how it depends on the operating conditions of the experiment. For clarity we 
specialize the discussion to to the Los Alamos proposal for a neutron edm search 
using Ultra-cold neutrons (UCN) and He 3 atoms diffusing in superfluid He 4 as 
a co-magnetometer, ^2] but the generalization to other cases is straightforward. 

First analyzed by Commins |5] in the context of a beam experiment, the 
frequency shift has been discussed in some detail by Pendlebury et al JU| in 
connection with experiments involving stored particle gases. Additional discus- 
sion and calculations have been given by 

Our present understanding of the effect can best be summarized by figure 1, 
which appeared as figure 3 in . This is a plot of the normalized (linear in E) 
frequency shift vs. normalized Larmor frequency for various values of collision 
mean free path and wall specularity. 

These results have been obtained by numerical simulation of the position- 
velocity correlation function and taking the Fourier transform. According to 
[[TT]. equ 26] the frequency shift is given by 



where R(t) is the position velocity correlation function defined in equ. 



appealing to try to make use of the zero crossing, apparent in figure 1, to reduce 
the effect. 

In this note we present an analytic model for the correlation function R (r) 
and compare it to the results obtained previously by numerical simulations. 
Using the limiting form of this model, valid for collision mean free paths small 
compared to the vessel size, we calculate the temperature dependence of the 
frequency shift for 3 He diffusing in superfluid He. This is an experimentally 
favored regime as collisions are seen to reduce the effect, as well as the slope at 
the zero crossing, considerably. We also propose a method for measuring the 
spectrum of the correlation function, i.e. the frequency dependence of the shift, 
directly. 



2 Analytical model for the correlation function 



2.1 Gas collisions 

1) We consider a particle that moves among scattering centers. If r c is the 
average time between collisions the velocity autocorrelation function will have 




(1) 



27]: 



and a = 
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the form [H| 

ip(t) = (v(t)v(0)) = v 2 e~^. (3) 
In the other words, ij)(t) obeys the equation: 

a (4) 

at t c 

2.2 Wall collisions 

We consider particles moving in a cylindrical storage cell in a vacuum. As 
shown in |10 | .|11 | the frequency shift depends only on the motion in the x, y 
plane. Referring to figure 2, the trajectory sweeps out an angle 

a = arccos (r/R) 

with respect to the center in a time 

T wa ii \/R 2 -r 2 



where T wa u is the time between wall collisions. The average angular velocity for 
a single trajectory is then 

arccos (r/R) v 

uj = = — . 

VR 2 - r 2 

The average squared frequency for all particles with velocity v (in the x, y plane) 
has been given by ^UJi e q u i (28) as 



As the motion is not strictly circular each trajectory will experience a compli- 
cated spectrum for the time varying field (Ref. JU|> section IV D and figure 7), 
but the mean square of the fundamental frequency will be given by © • 

Assuming that the result is dominated by the fundamental frequency, the 
velocity autocorrelation function will, to this approximation , obey the equation: 

iw + (j2 )m=0t (6) 



2.2.1 Non-specular wall collisions 

Following the proceeding section we are modelling the correlation function as a 
harmonic oscillator. During one traversal of the cell the oscillator will undergo 
a phase change 

<p = uiT W aii = 2a = 2 arccos (r/R) 
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A non-specular reflection from the wall would result in a change in the incident 
angle for the next collision, x, by a random amount A% and hence a change in 
the accumulated oscillator phase by 

= ?A X 
dx 

d<j) d<t> dr 
dx dr dx 



With 



sinx — r/R 



we have 



cos x = yl - r 2 /R 2 

d X = 1 

dr Vi? 2 - r 2 

d<f) 2 



so 



dr Ry/l-r 2 /R 2 



A0 = 2A X 



Since the changes A0 are random the phase <fi will make a random walk so that 
after a time t we will have 

(A0) 2 \ =4((A X ) 2 \^ = 2((A X ) 2X "• 



Twall \ I \JR 2 - r 2 

Averaging the amplitude of the oscillator over the distribution of A<fi the am- 
plitude will be reduced by 

(cos0) = l-i((A0) 2 )^l-— exp / 



t 7~non—spec \ ' non-spec 

Thus we have 

1 



(Ax) 2 )« 



Averaging over r we find 

^\ = <(A x) >/*Z^ = i((A x >>}» 

Tnon-spec / \ ' JO V R 2 - T 2 7T \ /it 

where i* 1 (r) dr is the probability that a trajectory has a pericentric distance 
between r,r + dr : 

. 4Vi? 2 - r 2 
F(r) ^^i?^- 
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2.3 Model expression for the correlation function 

When the particles are moving inside a vessel containing scatterers, one can 
expect that the equation for the velocity correlation function is the combination 
of J3J and ©, thus we are led to write 

This is just the equation for a damped harmonic oscillator and its general solu- 
tion is of the form: 

V>(t) = cie-" lT + C2e-'» T , (8) 
where 



Then, we have for the boundary condition satisfied by t/)(t), 

T 

h(r) = J ip(t)dt = — ^— — ► 0, when t — * oo. 
o 

Thus, the velocity correlation function [ip(r) — (v(t)v(Q))] will have the form: 
M T ) = mv2 ( e -mr _ Hl e -V2r\ ( 1Q j 

771 — 772 V m ) 

and the function R(t) will be given by: 

r ^ 

R(t) = 2 1 Mx)dx = — — f 1 - e'^-^A e -" 2T . (11) 

7 77l — 779 V / 



2.3.1 Overdamped, short mean free path, limit 

In the overdamped limit 



we find: 



Therefore: 



27 » ^0 (12) 



- — J ^2 ^ T c (o>q) , 77! > 77 2 . (13) 



R(t) = 2Xv (1 -e~£) e"*, (14) 



with 



using (JSJ). 



T -^>-(^)l\- (15) 
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The numerical simulations of the correlation function shown in figure 1 of 
[TT] (for a cylinder) , were seen to be in good agreement (in the short mean free 
path limit) with equ. $Tty above (equation 43, ^J) with 

T = .6R 2 /Xv (16) 

where A is the collision mean free path, v is the particle velocity and R is the 
radius of the cylindrical measurement cell, (taken as R — 25, in the rest of this 
paper). The factor 0.6 in equation (|16|l is in good agreement with the factor 
(6/tt 2 = 0.608) in equation (jT5|) . 

In the case of He 3 moving in superfiuid H e 4 , t c , v and A are all functions of 
temperature For UCN on the other hand the only collisions are those with the 
walls, all parameters are independent of temperature and we have to take the 
limit 

In figure 3 we show a comparison of the correlation functions calculated by 
numerical simulation, red, the model equ. (| 1 II) . blue and the short mean free 
path limit 1|14|) . green. 

In the figures r a = R/X. We see the model gives a reasonable fit for a wide 
range of r and the short mean free path expression (|14J) . (green), is a reasonable 
fit for r > 2.5. The disagreement for smaller r is most likely due to the neglect 
in our model of higher harmonics of the motional frequency 

In the following we will use (|14|l to calculate the frequency shift averaged 
over a Maxwell-Boltzman velocity distribution as a function of temperature for 
3 He diffusing in 4 He, taking into account the velocity dependence of the mean 
free path. 



2.4 Calculation of the zero crossing frequency for the vxE 
induced systematic 

According to JQl the frequency shift is proportional to S (oj), the cosine Fourier 
transform of R(t) which is given by J2J). For exp(— a|r|) the cosine Fourier 
transform is a/ (a 2 + oj 2 ) so that we have, using the short mean free path limit 
of the model, equation l|14|) 



S(oj) 



Vl + V2 



+ V2 OJ 2 + (f]i + f]2Y 



Xv 



(17) 



This is easily seen to have the correct limits for u> — > 0,oo, (equations 54 and 
70 in 21]). To find the value of oj where the effect goes to zero we put 

S (oj) = 

OJ + 771 + 772 = (w +%) 

T)2 
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in the case where r]i ^772 (A -C R) . Thus we see that the zero crossing, in 
the limit, is at 

/ Aw v 

which should hold for reasonably small values of X/R. Thus for short A all 
curves should have the same zero crossing. This is seen to be satisfied by the 
numerical simulations of figure 1. 

3 Frequency shift averaged over velocity distri- 
bution 

Using the limiting analytical model (|14f> for the correlation function and its 
Fourier transform we can perform an average over the velocity distribution. We 
do this for the case (1) of a constant collision man free path and then for the 
realistic case (2) of a mean free path proportional to velocity in order to show 
the influence of a variable collision mean free path. 

3.1 Mean free path independent of velocity 

The frequency shift, is given in terms of the spectrum S(u>) l|17l) by: 

du)(u>,T) = abS(uj) 



Suj(uj,T) = abR 2 (.6) 



(fc(T)/.6) 2 



(fc(T)/.6)(fc(T)/.6+l/fc(T)) 



(ivR/v) 2 + (fc(T)/.6) 2 (u;R/vf + (fc(T)/.6 + l/fc(T)) 2 

(18) 

with k = X/R = l/r a being a function of temperature (A(T) = 3D(T)/v(T)) as 
shown in figure 3c), (D(T) = 1.6/T 7 is the diffusion coefficient for He 3 in He 4 
as measured by and we have taken R = 25) . 

The average of this over the velocity distribution of the He 3 will be given 

by 



B 2 



= abR 2 (.6) 



1 



a 2 /ai 



(u>R/a 2 ) +v 2 (uR/a^+v 2 
where a 2 = (fc(T)/.6), cti = (fc(T)/.6 + l/fc(T)) . Now we have 



v A e~ v 1$ 



dv 



a 2 + x- 



: e- x2 ^ 2 dx ee F(a, (3) = 2a 3 V^e^ f 1 - erf 



- 2a 2 (3 + (i 3 
(19) 
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so that 



(M^)> = — (-6) 



abR 2 
abR 2 n(u,T) 



i ujR , A eta ( ujR , . 
F ( — ) F I — ,P(T) 



where /3(T) = 1.28 x 10 4 \j j^cm/ sec is the most probable velocity for He in 

He 4 . F{a,(3) tends to be difficult to calculate numerically for a/ (3 > 5 so we 
will use the asymptotic expansion (Ref. section 7.1.23, pg 298) 



fire z erf c{z) 



z { 2z 2 + 4z 4 8z 6 



for a/j3 > 4. Note that the leading term in the expansion is canceled by the 
term 2a 2 /3 in (|19|l . To plot the results we choose 

X (lu,T) = luR/(3{T) (20) 

as the independent variable. 

Figure 4 shows the normalized, averaged over the velocity distribution, fre- 
quency shift plotted versus the frequency relative to the appropriate, most prob- 
able velocity, pp")l: 

Expanding the plot we see the region near the zero crossing (figure 5). 
and expanding still further (figure 6) 

In this figure we see the zero crossings converging as the mean free path gets 
smaller (increasing temperature) just as in the single velocity case fsection [2.4J) . 

From the point of view of an experiment it is perhaps more interesting to 
plot the results as a function of temperature. 

Figure 7 shows the temperature dependence of the velocity averaged fre- 
quency shift for a range of frequencies specified by the corresponding values of 
X (u, T), (equ. |23, i.e. u = ^X (w, T). 

Expanding the plot (figure 8) shows the region where the effect can be min- 
imized: 

Here we see that the applied frequency can be chosen so that the velocity 
averaging does indeed reduce the effect. In fact reductions of 10 -5 to 10~ 6 
appear feasible. 



3.2 Mean free path oc velocity, cross section ~ 1/v 

We have k(T) = X(T)/R. Since the velocity of the He 3 is much less than the 
phonon velocity (2.2 x 10 4 cm/ sec) the collision rate of phonons with the He 3 
will be independent of the He 3 velocity. In a time r c a He 3 with velocity v, will 
move a distance X V (T) = vt c (T). Thus 

k v (T) = X V (T)/R = vt c {T)/R = v/s{T) 
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with s(T) = R/t c {T). 
From equation (|18|) 



5u(u},T) =abR 2 (. 



(k v (T)/.6) 2 



(fc u (T)/.6) {k v (T)/.6 + l/k v (T)) 

(uR/v) 2 + (k v {T)/.6) 2 {cuR/vf + (k v (T)/.6 + l/k v (T)f 

(21) 



The average of this over the velocity distribution of the He will be given by 



(Su(u,T)) 



dv 



= abR 2 (.6) 



abR 2 ^ (w,T) 



{cjR/a 2 f + v 2 (uR/aif + v 2 



/? 3 



(22) 
(23) 



where «2 = (k v (T)/ .6), otx = (k v (T)/.6 + l/k v (T)) . In the region of interest for 
suppression of the effect, k v (T) <C 1, we have oc\ ~ l/k v (T). Then the integral 
in (|22ll can be written 



V 



1/.6 



( -6s V / ^fl 



where y = v/(3(T). The second integral can be evaluated as 



4 " y dy 



y e 



Defining 
and writing 



4 — y J 3 

we y dy = 

2 4 



s(T) 



R 



p(T) 
R 



0(T) P(T)t c (T) A c (T) fe(T) 

where A C (T) and fc(T) are evaluated at the most probable velocity j3(T), we 
have a(T) = (^) (^) = ^yx(^T) 

The first integral can be evaluated in terms of a hypergeometric function 
(hypergeom ([1] , [jj — j] , ~j x2 )) but the series for this has some convergence 
difficulties for the parameters of interest. Thus we define 



/(«) = 



■e y dy 
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(to be evaluated numerically as the integrand is well behaved) so that 



*(w,T) = (.6) 



k(T) 



x(u,T) 



2.5 



(x( W ,T)) 2 +(^) 



(24) 



In figure 9 we plot the frequency shift averaged over velocity as a function of 
Larmor frequency for fixed temperature taking into account the velocity depen- 
dence of the mean free path. 

Figures 10 and 11 show the frequency shift as a function of temperature for 
various frequencies 



4 Measurement of the correlation function, R (r) 

While the theory presented here is expected to be accurate for the case of 
polarized He 3 diffusing in superfluid He 4 it would be nice to be able to confirm 
it experimentally and to check its applicability to other edm searches in progress. 
In this section we present a generally applicable, straightforward method to 
measure the frequency spectrum of the relevant correlation function (|2J) and 
hence the frequency dependence of the frequency shift directly. 

If we apply a uniform magnetic field gradient dB z /dz large enough so that 
it dominates all other field gradients that may be present, there will be a radial 
field in the x - y plane S r = - (~r /2) (dB z /dz) . Then 

dB x , y ldB z 



dx, y 2 dz 
and the field correlation functions will be 

(B Xi (t) B Xi (t + t)) = (af (an (t) x, (t + t)) 

where Xi — x or y. Then following McGregor |14| the relaxation time, Ti will 
be given by (in the case when dB z /dz is large enough so that wall relaxation 
can be neglected) 

with 

/>oo 

S r (oj)= / (¥x (t) ■ r± (t + t)) e- lUT dr 



/OO 
(r) coscurdr 
-oo 
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where we recognize that the correlation function Rj?j* (r) is an even function 
of r. A measurement of T± will thus yield the function S r (uj). Now 



where we used (r) = — cP-R-^y (r) /c?t 2 ^3 f° r the velocity correlation 

function. Then 



By comparison with equ (38) of we identify u> 2 S r (oj) /2ir with ip (oj) of that 
work and we then see that the systematic frequency shift is given by (equation 



Since we can determine S r (oj) from the T\ measurements we can determine the 
frequency dependence of the systematic frequency shift for any experimental 
conditions without the need of applying an electric field. 

5 Arbitrary Magnetic field geometry 

Our discussion has assumed a magnetic field configuration with G z — dB z /dz 
constant. Pendlebury et al jlllj have shown, using a geometric phase argument, 
that regardless of the field geometry the effect only depends on the volume 
average of G z in the high frequency (called by them the adiabatic) limit. In a 
recent note, Harris and Pendlebury ^7] have shown that in the case of a field 
produced by a dipole external to the measurement cell, this does not hold in the 
low frequency (diffusion) limit. In this section we discuss this problem using our 
correlation function approach in order to give some physical insight into what 
is happening and display details of the transition from one case to another. 

5.1 Short time (high frequency, adiabatic) limit of the cor- 
relation function 

Reference [TI], has shown that the systematic edm is given, in general, as the 
Fourier transform of a certain correlation function of the time varying field seen 
by the neutrons as they move through the apparatus. Equation (23) of that pa- 
per gives the frequency shift proportional to E as ( oj (t) lies in the x, y plane) 





40 of mi) 




(25) 



6u> E (t) 



1/ dr 

2 Jo 



{ 



COS UJ a T [tO (t) X ~UJ (t — t)] 

+ sinw T [u> x (t) lo x (t - t) + Lo y (t - r) u y (t)} 




12 



It can be shown that the term multiplying sin uj q t goes to zero on averaging 



over a uniform velocity distribution ((v x v y ) = 0, 
V x B — O.Then, for short times, r, 



v /2J and using 



Stu(t) 









C0SUJ o T 












COS UJ T 







, _> , , d uj Id 2 u) 



/dTo 



1 d 2 ~uj 



2 dT 2 



(27) 



We are considering values of r so small that the velocity doesn't change in 
that time interval (r < t co «) . 



Then 



d~uj 
~~dt 
d 2 U 
~d^ 



7 (~B xy {t) + lT/c x £f) 

7 (v5 (1? (*)) • V 
^ d 2 v3 
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and 
Suj(t) = 



1_ 

2c 



dr cos uj t 



- (l3 xy {t) +1? x i?) X ( 



ld 2 z3 



dTo 

~dT T ~2~ch 2 



(28) 



The term linear in E and r is then 



5a; (t) = — — / dr cos uj t 



Vi? • » rX t) X E 



defining 



1 E f* 

—— / dr cos uj t (ar) 
2 Jo 



(29) 
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We have now calculated the correlation function for short times. It starts at zero 
at t = and rises as ar. Eventually it will reach a maximum. By concentrating 
on the high frequency (uj ) behavior of Suj the result will be independent of the 
details of the maximum, depending only on a. Thus we can replace ar in l)29|l 
by sin ar or any function with the same initial slope. Thus we are led to take 



Suj(t) = lim / dr 

2 w„^ooj 



lim 



~/E 



2 uj 2 — a 



cos uj t sin ar 
jE a 



E 
~2c~Eft 
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Introducing components taking averages and using V • B = this reduces to 

1 I dB 



in agreement with equ. (2) of ^5 lr i m that equation, R 2 ui 2 is replaced by 
(w 2 ) = v 2 /2. We have shown that in the adiabatic (short time) limit the 
systematic (false) edm effect depends only on (^7-) regardless of the geometry 
of the magnetic field, a result obtained previously by Pendlebury et al ^U] and 
confirmed in |17| . 

The next order term in (|28|) is easily seen to be of order v 3 t 2 and so will 
average to zero, the next term which contributes will be of order v t 3 and so 
will be negligible in the short time limit we are considering. The condition for 
this to be valid is (vt/A) 2 <C 1 where A is the scale of variations in the applied 

magnetic field (jfe-j^ ~ A" 1 ) 

5.2 Longer time behavior of the correlation function 

For long times the expansion 127|) is clearly not valid and we must expand in a 
series in the spatial coordinates. We start from 



7 2 

8uj - 



J drcosuj T (B' (t) xB'(t- r) 



2 

where b = E/c, the brackets represent an ensemble average and 

B' x = B x (?(t))-bv y 
B' y = B y (*)) + bv x 

Then we write 

(3' (t) x ?> (t - r)) = b ( % %(t)) «b (* - rY - v y (t) B y (? {t - r)) 
\ Vj \ — [B x (r [t — t))v x (t) —v y (t — t) By (r (t))) 

(31) 

and expand the field in a Taylor series 



B x (V (t)) 



B x (0,0, 0,t) 
d 2 B 



dB x 



Ox 



x(t) 



8B X 



dy 



dx 2 

)2 , 



f d 2 B x 
\ dxdy 
fd 3 B x 



V dx 3 



x 2 (t) + 
y{t)x(t) 
x 3 (t) + 



d 2 B x 



dy 2 
, d 2 B x 



y 2 (t) 



y(t) 

o 

d 2 B x 



dB x 



dydz 
d 3 B, 



dz 2 
z (t) y (t) + 



dz 



z 2 (t) 



d 2 B,, 



dzdx 



••(*) 



X (t) z (t) 



dy 3 



y 6 (t) 



14 



(similarly for B y ). Concentrating on the first and last terms in (|31l) and noting 
that there are no correlations between any functions f(xi,Vi) and g (xj,Vj) we 
find 



E 

Xi—x,y 



Xi (t) 



dx] 



A (t) 



d 3 B x 



dx 3 



4 {t)...)v Xi (t-T)-{(t)^(t-T)} 



where the second term is obtained from the first by interchanging (t) and (t — r) . 

By symmetry we see that (x\ (t) v Xi (t — r)) = so that the next contribut- 
ing term will be proportional to 



d 3 B :J 



dx~ 



(xf(t)v Xi (t-r)> 



The first order term will be proportional to 

dB x , dB v 



We see that the condition 

d 3 B Xi 



dx 3 



dx 



R 2 < 



dy 
dB x% 



dxi 



dB z 
dz 



R2 1 
or — «1 



will insure that the higher order terms can be neglected. In the extreme case 
considered by Harris and Pendlbury, |17| . this condition is strongly violated so 
our method cannot be applied since the higher order terms remain significant. 



6 Discussion 

We see that in both the case of a velocity independent mean free path and 
one proportional to velocity the behavior of the systematic frequency shift is 
essentially unchanged. What is changed is the frequency where the effect goes 
to zero, i.e reduced frequency to' ~ 1.6 in the case (1) of a constant mean free 
path and cu' ~ 2 for the case (2) of a mean free path proportional to velocity. 
This is because longer mean free paths contribute more to the frequency shift 
as seen from figure 1, and thus higher velocities contribute more to case (2). 
This shifts the effective frequency to higher values for the same value of reduced 
frequency. 

Due to the heavy mass and slow velocity of the He 3 , Baym and Ebner [T5| 
conclude that the phonon scattering on He 3 is predominantly elastic. Single 
phonon absorption is kinematically forbidden on a single He 3 and can only 
take place as a result of He 3 — He 3 collisions which will be negligible for the 
low He 3 densities considered here. Thus our approach, where we calculate the 
correlation function for an ensemble of trajectories with constant He 3 velocity 
and then average the result over the velocity distribution should be an excellent 
approximation. 
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Note that for temperatures T > .38if, r Q = 1/k > 20, so that according to 
fig. 3 the limiting expression for the correlation function should be quite accurate 
and figures 10 and 11, which are based on an exact average over the Maxwell 
distribution of the He 3 velocities, imply that one should be able to control the 
effect to high degree. 

What would be ideal, would be to set the size of the shift for He 3 to be 
slightly negative i.e. equal in sign and magnitude with the shift for UCN. The 
UCN, with their relatively low value of uj r {uj /uj r 3> 1), have a frequency shift 
which is represented by a point far to the right in figure 1). 

The limits on the ability to do this will be the usual experimental ones of 
the temperature and magnetic field gradient stability. 

We have shown that by varying the gradients to larger values it will be pos- 
sible to measure the spectrum of the velocity correlation function directly, thus 
allowing a precise determination of the frequency dependence of the frequency 
shift under the exact experimental conditions. 
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7 Figure Captions 

Figure 1. Note the curves are for a single fixed velocity. The velocity dependence 
is contained in the normalization of the frequency scale, uj r = v / R. 

Figure 2. Trajectory of a particle confined in a cylindrical vessel. 

Figure 3a. Comparison of correlation function, R (t) from numerical simula- 
tions and from the analytic model, red - numerical simulation, blue - complete 
model equation l|llfl . green - short mean free path limit equation (|14H where 
appropriate. r Q = R/X = 50, 10. 

Figure 3b) As figure 3a) r a = 0.5, 1, 2.5, 6, 7.5 

Figure 3c. Dependence of k = \(T)/R on temperature. 

Figure 4. Normalized, velocity averaged fequency shift vs. normalized fre- 
quency: X(uj,T) = ujR/P{T). red T=0.3K, blue 0.35K, green 0.4K, purple 
0,45K 

Figure 5. Normalized, velocity averaged fequency shift vs. normalized fre- 
quency, expanded scale, red T=0.3K, blue 0.35K, green 0.4K, purple 0,45K 
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Figure 6. Normalized, velocity averaged fequency shift vs. normalized fre- 
quency. Same as figure 8, but expanded further, red T=0.3K, blue 0.35K, green 
0.4K, purple 0,45K 

Figure 7. Velocity averaged frequency shift as a function of temperature, K, 
for various frequencies., red X=1.4, blue X=1.5, green X=1.7, purple X=1.8 

Figure 8. Velocity averaged frequency shift as a function of temperature, 
K, for various frequencies. As figure 7, scale expanded to show experimetnally 
interesting region, red X=1.4, blue X=1.5, green X=1.6, purple X=1.7, light 
blue X=1.8 

Figure 9. Frequency shift averaged over velocity distribution, allowing for 
velocity dependence of the He 3 mean free path, for various temperatures. X 
is the normalized frequency, X — uiR//3(T) with (3{T) the most probable He 3 
velocity, green T=0.35K, purple T=0.38K, red T=0.4K, blue T=0.43K. 

Figure 10. Frequency shift averaged over velocity distribution, allowing for 
velocity dependence of the He 3, mean free path, versus temperature, K, for 
various frequencies specified by the value of X — uiR/ /3 with j3 the most probable 
He 3 velocity, blue X=3, green X=2.2, red X=2, purple X=1.8. 

Figure 11. Expanded plot of the frequency shift averaged over velocity dis- 
tribution, allowing for velocity dependence of the He 3 mean free path, versus 
temperature, K, for various frequencies specified by the value of X — ujR/ (3{T) 
with (3(T) the most probable He 3 velocity, blue X=3, green X=2.2, red X=2, 
purple X=1.8. 
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